A transport coefficient: the electrical conductivity 
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Abstract. I describe the lattice determination of the electrical conductivity of the quark gluon 
plasma 0. Since this is the first extraction of a transport coefficient with a degree of control 

^— ^ . over errors, I next use this to make estimates of other transport related quantities using simple 

' kinetic theory formulae. The resulting estimates are applied to fluctuations, ultra-soft photon 

spectra and the viscosity. Dimming of ultra-soft photons is exponential in the mean free path, 

H ' an d hence is a very sensitive probe of transport. 
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ly-) , 1. Introduction 

Transport coefficients have recently been computed to leading logarithmic order in weak coupling 
perturbation theory ^ Ej. However, at temeperatures of interest to heavy-ion collisions 
Qh! experiments, the gauge coupling g = 0(1) and this computation is ineffective. At the other 

extreme of g — > oo, the AdS/CFT correspondence has been exploited to bound the shear viscosity 
rj/S > 1/4-7T [3j, where S is the extropy density. This result is for a highly supersymmetric version 
'"^ \ of QCD and needs to be corrected in order to connect with experiments. 

There is clearly a case for a lattice computation, if one can have control over the process 
of analytic continuation from Euclidean to Minkowski space-time. Direct computations of the 
shear viscosity were tried (HE], but were inconclusive because statistical problems beset even 
the Euclidean computations of the correlator of the energy-momentum tensor, thus making it 
hard to check the analytic continuation. 

In this work |B] I shift the focus slightly and compute the electrical conductivity, a, of the 
quark-gluon plasma. This involves a lattice computation of the correlator of the electromagnetic 
current. Since this is no harder than the extraction of the p-meson mass, statistical errors in 
the Euclidean computation are under control, and effort can be concentrated on the harder job 
of analytical continuation. We are able to demonstrate control over this process, and thereby 
extract a. 

This is interesting in itself for heavy-ion collisions because several properties of soft photons 
in media can be inferred once a is known. Since this is the first statistically reliable extraction of 
a transport coefficient from lattice QCD, one can use this to estimate many different quantities 
which depend on transport, and are important for heavy-ion physics. With a in hand, and using 
two well-specified approximations, we are able to predict the value of rj, and compare the lattice 
results with others. 

The lattice computations described in Section 2 are from The material in Section 3 is 
partly new. 



2. The lattice computation 

All lattice computations of transport coefficients are performed in the context of linear response 
theory. The response, A, of a system to a time-dependent force F is given, in the frequency 
domain by A{uj) = x{u)F{oj). The response function x( w ) can then be written in terms of the 
retarded correlation function of A, and the transport coefficient can be extracted using Kubo 
formulae. For a, the Kubo formula is 

(1) 

uj=0 

where the sum is over spatial polarisations of the imaginary part of the retarded photon 
propagator, i.e., the spectral density, p E v M , of the electromagnetic current correlator. The 
usual QED Ward identity gives p 00 (uj, 0, T) = 2ttxqoj5{uj) = 0, where \Q is the usual charge 
susceptibility. 

A lattice computation proceeds from the spectral representation for Euclidean current 
correlators — 

roc A, j 

G EM (t,T) = —K(Lj,t,T)p BM (u,T), (2) 

JO 27T 

where K is the free propagator, G EM is the product of the zero (spatial) momentum vector 
correlator summed over all spatial polarisations, G v , and the EM vertex factor C EM = 
AnaJ2f e 2 f, where e/ is the charge of a quark of flavour /. For later use we note that C EM 1/20 
for two flavours. 

The main problem in this determination is that p EM is needed at an infinite number of points 
and G EM is known only at Nt different values. The solution is to constrain the function p EM 
through an informed guess, and use a Bayesian method to extract it. 

We note a second complication. A free field theory computation shows that p EM (u,T) 
increases as u> 3 at large co. This means that a naive application of the Kubo formulas would 
violate causality We try to solve this problem by using AG EM (t,T) which is the difference 
of the measured and free-field theory propagators to extract Ap EM . 



*(T) = ±—p EM i(u,0,T) 




Figure 1. The spectral density Ap EM for 
T = 2T C on a 12 x 26 2 x 48 lattice with quark 
mass 0.03T C . See the text for a discussion of 
statistical and systematic errors. 



Figure 2. The propagator AG EM obtained 
through the fit of Ap EM compared to the data 
for T = 2T C and Nt = 8 (circles), 10 (squares) 
and 12 (pentagons) with quark mass 0.03T C . 



Using a variety of Bayesian methods [0] to extract AG EM (t, T) we found a peak at io/(2nT) ~ 
0.15 (Hj and a smooth drop to zero at u = (see Figure^). The checks which have been performed 
on this extraction are — 



(i) The statistical errors are negligible, as seen in Figures ^ and [2J 

(ii) The position of the peak seen in Figure ^ is not strongly dependent on the lattice spacing, 
as checked by using Nt = 8, 10 and 12. 

(iii) Variations in the results are investigated due to algorithmic changes in the Bayesian prior 
and the truncation and discreitzation of eq. ©• The width of the band of systematic errors 
shown in Figure ^ is thrice the width due to such variation. 

It would be useful to have similiar tests of robustness on all uses of Bayesian methods including 
the maximum entropy method (MEM). 

The goodness of fit to the Euclidean data can be judged from Figure [21 We draw attention 
to the small error bars on this data when compared to the errors on typical measurements of 
Euclidean correlators of energy momentum tensor. 

Having established the existence of the peak, further information is extracted using a 
parametrization of the peak in Ap EM as the ratio of two polynomials, with that in the numerator 
being of lower order than the denonimator, so that the integrals used in the Kramers-Kronig 
relation converge. The denominator has to be at least a quartic polynomial in order to smoothly 
connect to the weak-coupling picture of a pinch singularity giving rise to the peak along real 
us jHj- (We note in this connection that earlier attempts to extract transport quantities have 
worked with the full p EM rather than with Ap EM . This error can strongly suppress estimates of 
the transport coefficients. Since p EM increases with to, this biases the fitted peak towards larger 
u, thus decreasing the slope at the origin for any fit form with one peak.) Using these fits our 
estimate of the electrical conductivity is — 





f 7.5 ±0.8 


(T = 


1.5T C 


C E M X | 


7.7 ±0.6 


(T = 


2T C ) 




[ 7.0 ±0.4 


(T = 


3T C ) 



3. Appplications 

3.1. Viscosity 

In kinetic theory one can write the Drude formula for the electrical conductivity — 

a = CbmS « t \ (4) 
m 

where S q is the entropy density in quarks, m their screening mass and r q their mean free time. 
Using the lattice result for a one gets r q = 0.30 ± 0.03 fm. While charge is transported through 
quarks, momentum is transported mainly through gluons. Due to colour factors the mean free 
time for gluons is about half the mean free time for quarks. Then using a kinetic theory formula 
for the viscosity along with such an estimate of r 5 one gets 

\ « 0.21 ±0.02 > — . (5) 

S 47T 

Only lattice errors are shown in this estimate. Errors on the kinetic theory are harder to extract. 
The lattice + kinetic theory result is 2.5 times the AdS/CFT limit. 



3.2. Photon dimming 

Non-zero electrical conductivity leads to screening of photons with energy below the peak of the 
spectral function found earlier, i.e., for u < 300 MeV. This is a direct consequence of Maxwell's 
equations. This range of screened energies lies far below the excess which may have been seen 
by the CERN NA49 experiment. It seems possible that direct photons at such low-energy could 



be distinguished from photons coming from decays of pions at RHIC-II. One method for using 
an observation of photon dimming to extract the conductivity has been suggested 
The mean free path, £, of these ultra-soft photons is given by the formula 

I = « 6 fm. (6) 

Since the fireball at RHIC is about 7 fm in size, one could expect dimming by a factor of 
exp(— 7/6), i.e., by a factor of 3.2. There is a self-consistency argument here involving only the 
data. If the viscosity is small, then mean free times of quarks and gluons are small, and the 
mean free path of photons can only be about 20 (= 1/C EM ) times larger. However, dimming is 
exponential in the mean free path. For example, the mean free path as predicted in AdS/CFT 
would lead to dimming by a factor of 18.8. Thus, the use of photon dimming to extract a is an 
important cross check of values of r\ extracted by other means. 



3.3. Fluctuations 

The study of event-to-event fluctuations in limited rapidity bins is expected to give information 
on the phase diagram of QCD as well as make contact between experiments and basic lattice 
QCD computations. However, fluctuations of conserved quantities in a limited region of space 
are expected to even out by transport processes such as diffusion. The computation of a is 
relevant not only to the diffusion of charge, but also of baryon number and isospin. 
The diffusion constant can be obtained using the kinetic theory formula ^U] 

D = r q c 2 s 0.1 fm, (7) 

where we have taken the square of the speed of sound c 2 s = 0.30 ± 0.01 at T = 2T C , as shown 
by recent lattice computations If the plasma is highly dissipative, then D must be small. 
When D decreases, the fluctuation are more localized at the same time. Small viscosity implies 
smaller mean free path, and hence longer persistence of the fluctuation signal. In the approach 
of ^0] this would mean that the rapidity width of a fluctuation is AY < 2D/tq ~ 0.4 for 
fluctuations created at time tq = 0.5 fm. 



3.4- A brave new world 

When weak coupling was the only tool to investigate transport properties, it was felt that the 
QGP would be weakly dissipative. Now with other tools at our disposal, the question seems 
more open, and preliminary evidence is that dissipation is strong. The first lattice results on 
transport .6] pin down many aspects of transport, as outlined in this section, and lead to the 
hope that multiple experimental constraints on transport are feasible. 
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